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Abstract. In his seminal 1962 paper on the "threefold way" , Freeman Dyson classified the 
spaces of matrices that support the random matrix ensembles deemed relevant from the point 
of view of classical quantum mechanics. Recently, Heinzner, Huckleberry and Zirnbauer have 
obtained a similar classification based on less restrictive assumptions, thus taking care of 
the needs of modern mesoscopic physics. Their list is in one-to-one correspondence with 
the infinite families of Riemannian symmetric spaces as classified by Cartan. The present 
paper develops the corresponding random matrix theories, with a special emphasis on large 
deviation principles. 

Half a century ago, when physicists started to explore the usefulness of random matrix 
ensembles for the study of statistical properties of the spectra of heavy nuclei, their approach 
was firmly rooted in the classical framework of quantum mechanics. The Hamiltonian of a 
system was replaced by a random hermitian matrix each realization of which was supposed 
to commute with the appropriate unitary symmetries and with the correct "time reversals" , 
i.e. certain antiunitary operators. Concretely, the most general hermitian matrices that 
commute with time reversal in the literal sense are real symmetric matrices, whereas another 
type of "time reversal" symmetry leads to quaternion real matrices (see [121 Chapter 1] or 
[15] for details). In structural terms, the spaces of hermitian, real symmetric and quaternion 
real matrices can be viewed as tangent spaces to or infinitesimal versions of Riemannian 
symmetric spaces (RSS) of type A, AI, All, respectively. In his landmark article "The 
threefold way" ([9]) of 1962, Dyson proved that any hermitian matrix which commutes with 
a group of unitary and "time reversal" symmetries reduces to a block matrix, whose blocks 
are of the three types described above. 

In the 1990s the theoretical condensed matter physicists Altland and Zirnbauer argued that 
random matrix models for so-called mesoscopic normal-superconducting hybrid structures 
([I]; 0, [23] ) must be taken from the infinitesimal versions of the symmetric spaces of class 
B/D, Dili, C, CI (see the list of classical symmetric spaces in Section [1] below). What they 
had in mind are physical systems which typically consist of a normalconducting quantum 
dot, i.e., a small metallic device of extension lfim or less, which is in contact, via potential 
barriers, with two superconducting regions. The metallic quantum dot may or may not be 
disordered. In the latter case it is assumed that the geometric shape of the dot is such that 
the classical motion of a single electron inside it is chaotic. There may be a small magnetic 
flux present and some impurity atoms may cause spin-orbit scattering. The temperature is 
so low that the size of the quantum dot is much smaller than the phase coherence length 
of an electron. The latter feature is truly microscopic. On the other hand, tiny as it is, 
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the quantum dot consists of a huge number of atoms, which is a distinctively macroscopic 
feature. Systems of this kind, displaying both microscopic and macroscopic features, are 
called mesoscopic. 

We are not aware of an intuitive argument to explain why one should expect precisely these 
new "symmetry classes" (series of classical symmetric spaces) to arise in this context. But 
there is a more basic question we can try to give some hints about: Given that Dyson's 
theorem is a classification result seemingly in the most general setting quantum mechan- 
ics has to offer, why should there be any leeway for new symmetry classes in mesoscopic 
systems as above? The reason is that the standard account of the dynamics of this kind 
of system is phrased in the language of second quantization, i.e., linear superpositions of 
particle creation and annihilation operators are acted upon by Hamiltonians. For the meso- 
scopic systems described above, however, it is possible to convert this set-up into something 
akin to Dyson's first quantization framework (a Hamiltonian acting on a complex Hilbert 
space of "quasiparticle wavefunctions"). But there remains a crucial difference: The cre- 
ation and annihilation operators for fermions obey the canonical anticommutation relations 
c^cp + Cf3<J a = 5 a /3, c^ct + etc], = 0, c a C/3 + c$c a = 0, with a, /3 from an orthonormal basis 
of single-particle states, and these relations are mirrored in the translation as a symmetric 
bilinear form on the complex Hilbert space (along with the scalar product, of course). So 
one has from the outset more structural data than Dyson had, and it is no surprise that new 
symmetry classes arise. (A formalized version of this argument is in the first paragraph of 
Section [1] below. For details see the introduction to [T6]). 

Let us mention briefly that another source for new symmetry classes are Dirac fermions in a 
random gauge field background, leading to classes AIII, BDI, and CII of the list in Section 
□ below ([25], [T2], P31 Sec. 2.3 and 6.2]). 

After the physics examples for the new symmetry classes had been recognized, Zirnbauer and 
complex geometers Heinzner and Huckleberry took up the task of updating Dyson's classi- 
fication result to the new, enriched framework. In [16], they proved that the constituents 
which make up the Hamiltonian are in one-to-one correspondence with the ten infinite fami- 
lies of classical symmetric spaces as classified by Cartan (which motivates our choice of labels 
for the symmetry classes). 

Like Dyson's "threefold way" , the "tenfold way" of [IB] is established in geometrical terms, 
without reference to probability measures on the matrix spaces in question. It is the object of 
the present paper to provide a mathematical treatment of the corresponding random matrix 
theories. Besides being based on the systematic framework of [16], the present article differs 
from the existing literature on random matrix ensembles associated to symmetric spaces 
(see [7], [1]) in its focus on Large Deviations Principles. Thus, before going to business, 
let us review those aspects of the classical theory of the Wigner-Dyson ensembles which 
will be subsumed in the present analysis as instances of symmetry classes A, AI, AIL For 
convenience, we only mention the results in the GOE (AI) case. If X = {Xij)i<ij< n is a 
symmetric matrix of real-valued centred Gaussian random variables such that 

• (Ajj)i<j<j<„ are independent, 



then its distribution is a probability measure on the space of symmetric nxn matrices, which 
is invariant under conjugation by matrices from the orthogonal group O n . The eigenvalues 
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of X, Ai, . . . , A n , say, are real valued random variables, and by orthogonal invariance, their 
joint distribution has a Lebesgue density q n that can be given explicitly: 

1 1 n 

q n (x 1 ,...,x n ) = — II \ Xi _x il ex P(~2 n 5Z X ^ ^ 

n l<i<j<n i=l 

where Z n is for normalization. As n — > oo, the random measure 

1 n 
n *— ' 

tends to a nonrandom limit, namely, to Wigner's semicircle distribution with density 

il {W ^V2=F- (2) 

In their paper (3_J, Ben Arous and Guionnet have carried out a finer analysis of L n , and 
have proven that it satisfies a Large Deviation Principle in .Mi(R) (the space of probability 
measures on the Borel sets of M, endowed with the weak topology) with speed n 2 and good 
rate function 

7 ^ = \J x2 ^ dx ^ + \J J lo &\ x - y\~ x v( dx ) M*/) - ^ 

whose unique minimizer is the semicircle distribution. Recall that a family of probability 
measures (/i e ) £ >o on some topological space X is said to obey a Large Deviation Principle 
(LDP) with speed e~ l and good rate function I : X — > [0, oo] if 

• I is lower semi-continuous and has compact level sets Nl := {x G X : J(x) < L}, for 
every L e [0, oo[, 

• liminf £ ^ £ : log/U e (G') > — inf xe a VG C X open, 

• limsup e _ >0 elog/i e (y4) < — inf^g^ I(x) \fA C A closed. 

This paper is organized as follows: In Section [T] we review the symmetry classification of 
matrix Hamiltonians in mesoscopic physics, as given in [16J, and describe the infinitesimal 
versions of classical symmetric spaces that turn out to be in one-to-one correspondence 
with these symmetry classes. Then, in Section [21 we introduce probability measures on 
these spaces which enjoy invariance properties that reflect those of quantum mechanical 
observables. The resulting random matrix ensembles are called Hamiltonian ensembles. 
They can be viewed as generalizations of the Wigner-Dyson ensembles GOE, GUE, GSE. We 
use the geometric description of the underlying spaces to derive the induced joint eigenvalue 
densities in a uniform way, thus generalizing ([T]) above. In Section [31 we turn to the large 
deviations analysis of the empirical eigenvalue measure. We prove a generalization of the 
main result of [3], which covers not only the Hamiltonian ensembles, but also some matrix 
ensembles or particle systems of different origin that have been studied in recent years. 
In Section HI then, we describe the Gaussian Hamiltonian ensembles in concrete terms, 
analogously to the construction of the GOE which was reviewed above. We make explicit 
what the results of Section [3J mean in these special cases. 
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1. Symmetries of Nambu space 

Let (W, b) be a complex vector space of dimension 2n (n G N) together with a nondegenerate 
symmetric bilinear form b. By polarization we may assume that W = V © V* and that b 
is the natural pairing of V with its dual V*, i.e., b{x\ + <px,Xi + (^2) = ^1(^2) + ^2(^1)- 
Write S := f\V*. To if G V* assign the wedge multiplication operator e(<p) G End(S), 
and to v £ V assign the contraction operator l(v) G End(S). The e(cp) and l(v) satisfy 
the Canonical Anticommutation Relations (CAR) of creation and annihilation operators on 
fermionic Fock space, and if p : W — > End(S) is given by p{ip + v) — e(cp) + l(v), then 
(End(S),p) is a Clifford algebra Cliff(W, 6) for (W,b). We regard W as a subspace of the 
associative algebra Cliff (W, b), which we interpret in the usual way as a Lie algebra. One 
can embed so(W,b) as a Lie subalgebra consisting of elements which are quadratic in the 
w G W, and the adjoint action of so(W, b) on W turns out to be nothing else than the 
natural action of so(W, b) C End(PE) on W (see [13] for details). In physical terms, the 
embedding of W into Cliff(W, b) = End(S') suggests the interpretation of an element of W 
as a field operator on fermionic Fock space. The dynamics of a system of field operators is 
given by Heisenberg's equation of motion 

the self-adjoint operator H being the Hamiltonian of the system. Thus we have seen that 
if the dynamics of the system is governed by a quadratic Hamiltonian of a certain type, all 
relevant information is encoded in the structure (W,b). 

So far we have not yet made explicit that V, which plays the role of the space of single 
particle states, comes with an hermitian scalar product (•,•). It gives rise to a C-antilinear 
bijection C : V — > V* : v ^ (v, ■). C and (•, •) can be extended to the entire space W in 
such a way that {wi,w 2 ) = b(Cwi,W2) Vwi,w 2 G W. The triplet (W,b,C) is called Nambu 
space. 

Now suppose that a compact group Q acts on W by unitary or antiunitary transformations. 
It is inevitable to bring antiunitary transformations into play here, because time enters the 
formalism of quantum mechanics via ih4r, and so time reversal is an antiunitary rather than 
a unitary transformation. Once (W,b,C), Q and its action on W are fixed, the translation 
of Dyson's problem to the present context is as follows: Describe the space of all hermitian 
H G so(W,b) with the property that Hg = gH for all g G Q. It is convenient to call these 
H good Hamiltonians. Denote by Qq the group of all elements of Q which act by unitary 
transformations. The basic assumption of [16J is that Q is generated by Qq together with at 
most two elements acting by antiunitary transformations (One may think of systems which 
are invariant under both time reversal and charge conjugation). 

It is one of the key insights of Heinzner, Huckleberry and Zirnbauer in [TB], that it is possible 
to reduce the problem to the case Qq — 1. But this comes at the price that the good Hamil- 
tonians need not be elements of so(W, b), but can belong to 5p(VV, a) (for an alternating form 
a on W) or to sl(V) (diagonally embedded into End(V) ©End(V*)). Write s for any of these 
three Lie algebras. Although one is ultimately interested in hermitian operators, one first 
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considers the skew hermitian elements of s, which make up a compact real form g of s. Now 
let T be an antilinear transformation of W such that T 2 = ±id and suppose that Q = (T). 
From the reduction step in [16] it emerges that T may be assumed to fix g. If 9 denotes 
conjugation (in End(W / )) by T, then 9 restricts to an involutive Lie algebra automorphism 
of q. Let q = t © p be the decomposition into the +l-eigenspace t and the — 1-eigenspace 
p of 9, the so-called Cartan decomposition. Then ip consists precisely of those hermitian 
operators in q which commute with T, hence it is the space of good Hamiltonians. If G is a 
connected compact Lie group corresponding to q and K its closed subgroup with Lie algebra 
t, then p can be thought of as an infinitesimal version of the RSS G/K. If Q = 1, then the 
space of good Hamiltonians is iq. Since compact Lie groups can be given the structure of an 
RSS (see [15] . Ch. IV §6), this case fits into the overall picture. 



It is the main result of [16J that the following is the complete list of spaces of good Hamilto- 
nians that correspond to Nambu space with the kind of symmetries in question. The labels 
refer to Cartan's list of the classical compact Lie algebras and their involutive automor- 
phisms (see [TSl Ch. X §2.3]), i.e., plainly, to his list of classical symmetric spaces. Note that 
in order to obtain the case BDI below for the full range of parameters, and thus to obtain 
Cartan's full list from the symmetries of Nambu space, one has to refine the above argument 
in order to take care of the case Q = (T,Ti), 7\ being another antilinear transformation 
of W with Tf = ±id. We deviate from standard practice in Lie theory in that we do not 
require the matrices to be trace-free, in order to recover the familiar Wigner-Dyson classes 
as classes A, AI, AIL X* denotes the conjugate transpose of a complex matrix X. 
Class A: 

iQ = {X e C nxn : X hermitian} 



Class AI: 
Class All: 

Class AIII: 

Class B/D: 
Class BDI: 

Class Dili: 



ip 



ip = {X e R nxn : X symmetric} 

Xi X2 V Xt e C nxn ,X 1 hermitian, 
-X 2 X\ ) X 2 skew symmetric 



ig = {X E (iR) nxn : X skew symmetric} 



'PHI x* » : 



ip = <{ ( ^ ) : Xi £ (iR) nxn skew symmetric 



Class C: 



Xi X2 \ Xi G C nxn , Xi hermitian, 
X 2 — Xi J X 2 symmetric 
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Class CI: 

i P = | \ . X . e R »xn symmetric 

Class CII: 

where the space M. sxt of quaternionic matrices is embedded into C 2sx2i as 

M sxt = {( _y y J ■ U,VeC sxt Y 

More precisely, B/D splits into B for n odd, and D for n even. We will call classes A, AI, 
All Wigner-Dyson classes, BDI, AIII and CII chiral classes (in view of their role in modeling 
Dirac fermions) and the others superconductor or Bogolioubov-de Gennes (BdG) classes (in 
view of the discussion in the introduction, BdG being a keyword in the conversion to a first 
quantization set-up which was mentioned there). 



2. HAMILTONIAN ENSEMBLES 

In this section we randomize the good Hamiltonians, i.e., we put probability measures on 
Q := iQ resp. p := ip. Let G be a connected compact Lie group with Lie algebra g, K its 
closed subgroup with Lie algebra £. The adjoint representations Adc : G — > GL(g) and 
Ad^ : K — ► GL(p) are given by conjugation of matrices. If is a (nonempty open subset of 
a) finite dimensional Euclidian vector space, write m c for Lebesgue measure on o. 
Now we study probability measures on g resp. p, restricting our attention to those which 
are absolutely continuous w.r.t. nig resp. nip. Since quantum mechanical observables are 
invariant under unitary similarity transformations, it is natural to assume that the measures 
are invariant under conjugation with unitary matrices. So we only consider m -densities that 
are constant on the Ad^-orbits in q resp. nip-densities that are constant on the Ad#-orbits 
in p. 

Now we wish to compute the joint densities of the eigenvalues of a matrix leg resp. X G p 
chosen according to such a measure. These are easy consequences of the infinitesimal version 
of Weyl's integration formula for q and its analog, due to Harish-Chandra, for p. Here we use 
standard terminology of elementary Lie theory, see e.g. [22], [18], [8]. To state the formula 
for g, let T be a maximal torus of G with Lie algebra t, W = Na{T)/T the corresponding 
Weyl group, R + and t + C t compatible choices of a system of positive roots and of a positive 
Weyl chamber, respectively. 

Proposition 2.1. There exists c > such that for all f e L x (g, m ) which are constant on 
Ado-orbits there holds 

[fdm g = c [ f(T) H a{Tf m t+ (dT) 

f(T) J] HT)| 2 m t (dT). 




Proof. [8] Cor. 3.14.2 (ii) 



□ 
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To state the theorem for p, write a for a maximal abelian subspace of p, W = N^(a)/Ci^(a) 
for the Weyl group of g w.r.t. a, S + and a + C a for compatible choices of positive restricted 
roots (of w.r.t. a) and of a positive Weyl chamber, respectively. For p G £ + , m p denotes 
the multiplicity of p, i.e. the dimension of the joint eigenspace, corresponding to the linear 
form p, of the commuting symmetric operators which are induced on g by o. Alternatively, 
trip is the cardinality of the inverse image of p w.r.t the restriction process described in [18J, 
p. 263-4. 

Proposition 2.2. There exists c > such that for all f e L 1 (p,m p ) which are constant on 
Adx-orbits there holds 




Proof. PH Thm. 1.5.17] □ 

In practical terms, if / is an Ad^-invariant m^-density q — > [0, oo[ or an Ad^-invariant nip- 
density p — > [0, oo[, then f(X) only depends on X through its eigenvalues. Now, for all 
classes except A, AI, All, the nonzero eigenvalues come in pairs ±A. Note that in the chiral 
classes BDI, AIII and CII, the number of positive eigenvalues (without multiplicity) is s A t, 
and later on, when we will let matrix size tend to infinity, we will have to control the growth 
of s = s(n) as well as of n. What is more, note that many of the matrices that make up 
the spaces of good Hamiltonians are necessarily of even size, so that there is no question of 
simply letting matrix size tend to infinity. Hence we choose the following framework: n or 
(n, s(n)) is the (pair of) positive integer parameter(s) that will be controlled, t(n) := n — s(n) 
by definition. The space of good Hamiltonians for symmetry class C and parameters n or 
(n, s(n)), written as m£°, is contained in C d(n)xd(n) . The map 

7i n : M^ n) -> RPW.+, 

where 

R d '+ -={x eR d : xx > x 2 > • • • > x d }, 

assigns to each X e the nonincreasing vector (Ai, . . . , A p ( n )) of its (positive) eigenvalues 
(without multiplicity). By construction, ir n separates the adjoint orbits of G resp. K, and it 
is easily seen that 7t n (X) E RP( n )>+ a.s. 

Plugging well-known facts about (restricted) root systems (see, e.g., the appendix to the 
monograph [22]) into Propositions 12. 14 12.21 above, one obtains the following 

Corollary 2.3. Suppose that X is a random element ofM.^\ whose distribution is given 
by a Lebesgue density f such that f = f o -n n for some measurable f : ]RP( ra )' + — > [0, oo[. 
(a) If C is A, AI or All, then the joint density of the eigenvalues of x (in nonincreasing 
order) is 

(Xi, ... ,Xp( n )) >->• Const f(xt, ... ,X p (n)) ) f [Xi - Xj) P , 
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where (3 and p(n) are given in the table below. 
(b) Otherwise, the joint density of the positive eigenvalues of X (in nonincreasing order) 
is given by 

(xi, . . . , x p(n) ) i const /(ari, . . . , x p{n) ) \ | (x- - xj) /3 JJ s", 
where a, j3 and p(n) are given in the table below. 



Class 


d(n) 


p(n) 


a 




A 


n 


n 












2 


AI 


n 


n 












1 


All 


2n 


n 












4 


BDI 


n 


s(n) 


At 


n) 


\s(n) 


-t(n)\ 




1 


AIII 


n 


s(n) 


At 


n) 


2\s(n 


-t(n) 


+ 1 


2 


CII 


2n 


s(n) 


At 


n) 


A\s(n] 


-t(n) 


+ 3 


4 


B 


2n+l 


n 






2 






2 


D 


2n 


n 













2 


C 


2n 


n 






2 






2 


CI 


2n 


n 






1 






1 


Dili (n even) 


2n 


LfJ 






1 






4 


Dili (n odd) 


2n 


LfJ 






5 






4 



3. An LDP for the empirical eigenvalue measure 

In this section we prove a theorem which contains LDPs for the empirical eigenvalue measures 
of all Hamiltonian ensembles, assuming a product structure for / in the notation of Corollary 
12.31 Note that this subsumes the well-studied densities on of the form exp(— Tr(V(X))), 
V a polynomial with positive leading coefficient (see Remark 13.41 below) . This is the set-up: 
Let (3 > 0, 7 G N. Let S be a closed subinterval of R if 7 is odd, of [0, oof if 7 is even. Let 
(w n ) n( z^ be a family of continuous nonnegative real- valued functions on E. For 116N consider 
random variables A n = (A 1; . . . , \ p ( n )) with joint distribution Q n = P o A^ 1 G .Mi(£ p(n )) , 
given by its Lebesgue density 

p(n) 

q n (xi, . . . ,x p(n )) = — Yl \ x l ~ x]\ P Y[w n (x j ) n l j:p{ n)(x 1 ,...,x p{n) ), (3) 

n l<i<j<p(n) j =1 

where Z n is for normalization. In what follows we assume that p(n) — > 00 for n — > 00, 
satisfying 

lim EM = K e ]o,oo[. (4) 

Write J\f(f) for the set of zeros of a function /. We will make the following assumptions 
about the sequence (w n ) n : 

(al) there exists a continuous function w : S — » [0, oof such that 
- #J\f(w) < 00, J\f(w n ) C A/"(w) for large n. 
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— As n — > oo, w n — > w and \ogw n — > \ogw uniformly on compact sets, 

— logw is Lipschitz on compact sets away from Af(w). 
(a2) If £ is unbounded, then there exists n G N such that 

lim \ x \i«(fM)+e sup Wn ( x ) = o 

x— >±oo n>no 

for some fixed e > 0. 
For x — (xi, . . . , x p (n)) E £ p(n ) set 

. p( n ) 
L ^ := ^)XX 

Theorem 3.1. (P o (L n o A n y 1 ) n = (Q n o L~ l ) n satisfies an LDP on Mi(E) with respect 
to the weak topology with speed n 2 and good rate function 

m = ^ 2 [ /logi^-^rvwM^) (5) 



2 

— k I log w(x) jj,(dx) — c, (6) 



where \i e Ati(S) and 



c := lim \ogZ n < oo. (7) 



n^oo n 2 



COROLLARY 3.2. If I has a unique minimizer fi* , then 

F(L n oA n ^^*) = 1 

where — > means weak convergence. 

Proof of the corollary: Using the upper bound of the LDP one obtains the strong law apply- 
ing Borel-Cantelli's lemma, see [TQl Theorem II B.3]. □ 

Remark 3.3. If /3k > 1, then it follows from the theory of logarithmic potentials with 
external fields, applied to the weight function x i— >■ w(x 1//7 ), that // exists and that /(//*) +c 
is finite, see [21], Thm. 1.1.3 and Ex. 1.3.5. In this case (J7|) can be sharpened to |c| < oo. 
There exists a vast literature that describes \i* in more detail for various classes of weights, 
see e.g. [24], [5], [TH Thm. 3.1]. We will restrict ourselves to giving explicit formulae for /i* 
for Gaussian Hamiltonian ensembles below in Section HI In these cases, \x* is "universal" in 
the sense of Wigner's Theorem. This is the content of the companion paper [21]. 



Remark 3.4. Although the main focus of the present paper is on Hamiltonian ensembles 
associated to (infinitesimal) symmetric spaces, let us note that Theorem 13.11 contains an 
LDP for the much wider class of ensembles which is considered in [TTl (1.6)]. Furthermore, 
Theorem 13.11 applies to Wishart matrices (see Remark 14.21 below) and to Jacobi ensembles, 
for which the LDP was first proven in [20] . As observed in [7], the latter class of ensembles 
includes the random matrix ensembles associated to compact symmetric spaces, except those 
of classes A, AI, AIL 
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The proof of Theorem 13.11 extends the approach of [3], [H], [19]. Define onSxS the function 

F(x,y) := -^K 2 log|x 7 - y 7 | - ^(\ogw(x) + logtu(y)), 

with F(x,y) = oo for x 1 = y 1 (hence x = y by definition of E) or PI J\f(w) ^ 0, and 

its truncated versions 

F M (x, y) := F(x, y) A M, M > 0. 
Moreover, consider the functions 

F n (x,y) := -^(^i) 2 log |s 7 - y 7 \ - ^-^(logio n (x) + logtu n (y)) 

and their truncated versions 

F^(x,y):=F n (x,y)AM, M > 0. 
From the definition of _F n it follows that 

1 / 2n 2 n P(n) \ 

q n (x 1 ,...,x p(n) ) = —expl-—— F n {x i ,x,j) + ——^\ogw Tl {x l )\. (8) 

n \ yy ) i<i<j<p(n) p[ - ' »=i ' 

Lemma 3.5. 

(i) For any M > 0, F£f(x,y) converges to F M (x,y) uniformly as n — > oo. 

(ii) F zs bounded from below. 

Proof. The estimate log |x — y| < log(|x| + 1) + log(|y| + 1) implies 



F n (x,y)> 



2n 



log((|x 7 | + l)^ n)/n w n (x)) + log((|^| + l)M n)/n w n {y)) 



(9) 



Observe that log((|x 7 | + l)^")/* w n (x)j is bounded from above by (a2) and the continuity 
of w n . Invoking the full strength of (a2), one sees that for each M > there exist n 6 
N, -Rm > 0, S UjM > (v E Af(w)) such that F n (x,y) > M holds for all n > n on 

A M := {W V > R M } U [J {|x - i/| A |j/ - i/| < 8 V , M }. 

By compactness of A° M , (i) follows from (al) and the definition of F M ,F^. One also has 
F > Mon Am- F being continuous, it is bounded on A C M , and this proves (ii). □ 

Our strategy is to first consider the finite positive measures P n := Z n Q n . As to the upper 
bound, note that 

L n ®L n ({(x,y) E S 2 :x = y}) 



p{n) 

m RP („) -almost surely, since the eigenvalues are a.s. distinct under the product Lebesgue mea- 
sure. Hence almost surely 

J jf F^(x,y)L n (dx)L n (dy) = J J F™ (x,y) L n {dx) L n (dy) - -^j. (10) 
Now let A be a Borel set in A^i(E) and write 

A:= \xe £ p(n) : LJx) G a1. 
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Using (jHJ) and Holder's inequality we obtain 



l<i<j<p(ra) 

expl — — y^logw n (xi) )m RP{ n)(dx) (11) 

p(n) 

< J exp log w n (t)] m R (dt) 



(I) X (II). 



1/2 



Note that we have used that 2 Xli<j Fn( x i, x j) — Fn{ x ii x j) by symmetry of F n in its 
arguments. Now, lim^oo -4 log (I) = by (a2). On the other hand, for any M > (jTOl) 
yields 

(II) < j^exp (-2n 2 (l„(x)« 2 (F„ m ) - -^)) m MpW (^)' ' 



< <jexp( -2n 2 I inf /i 02 (F„ M )- ^ 



^ e A^ n p(n 
exp ( -n 2 inf ^ 2 (F n M ) ) exp ' 



1/2 



fj.eA J \V\ n ) 



Using Lemma [3. 5[ we obtain 

lim ( inf /i 02 (F n M )^l = inf p® 2 {F M ). 

n^ooKfi^A J n^A 

We have thus shown that for any Borel set A C A^i(S) one has 

limsup^logP n (L„ eA)<- inf ! ! F M \x,y) (i(dx)fi{dy). (12) 

rwoo n V&A J J 

Setting 

H(jj) := J Fdfi® 2 , H M (ii) := J F M d/i® 2 , 

one obtains well defined maps on A^i(E) (see Lemma I3T5T) . We will show that if is a good rate 
function that governs the LDP for (P n ). To this end, observe that, F M being bounded and 
continuous, H is weakly continuous on «Mi(£) for each M > 0. By monotone convergence, 
we have liniM^oo H M = H pointwise on A^i(E). As a limit of an increasing sequence of 
continuous functions, H is lower semi-continuous, i.e. the level sets {H < L} are closed. We 
claim that they are compact. Indeed, let mp '■= | inf F\ and a > 0. Then for any \x G .Mi(E) 
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one has 



inf (F + m F )(x,y) /i([-a,a] c ) < / (F + m F )(x,y) n(dx) n{dy) 

x,ye[-a,a] c J J J 

< H{n)+m F , 

hence {H < L} C K L , L e]0, oo[, with 



a>0 ^ V m a; 



-a,a]c(^ + ra F )(>,?/) 



Since lim^oo inf x ^_ a ^]c(F+m F )(x, y) = oo, is weakly relatively compact by Prohorov's 
theorem. Hence the rate function H is good. Furthermore, this argument easily yields the 
exponential tightness of (P„ o L^ 1 )^. Fix M > 0, and define K L (L > 0) as above, using F M 
in the place of F. For every \i G there exists a = a M > such that 

inf F M (x, y) + m F M ) //([—a, a] c ) 2 > L + m^a 
hence 

inf F^,y)/i([-a,a] c ) 2 >L. (13) 

x,y£[-a,a\ c 

Then (1121) implies 

lim sup \ log P n (L n EK C L )<- inf / P^/i 02 

M J.,®2 i / 77M J, ,02 



= - inf / F dfj, + / P M d// 

< - inf / inf P M (x,y) d/z 02 - inf / inf P M d/i 02 

< —L + m F M. 

Since inf P M > —00, we have shown that 

lim sup lim sup — logP n (L n £ K C L ) = — 00, 

L—>oo n—*oo Tl 

hence that (P n oI^)„ is exponentially tight. 

Now let -B(/i, 5) denote the ball centered at jj, G .Mi(£) with radius S for a distance compat- 
ible with the weak topology. Since /1 1— > H M ([i) is continuous, from (JT2l) we obtain for any 



inf lim sup -i log P n (L n G P(/i,5)) <~ f f F M i 



5 

Finally, letting M go to infinity, we obtain the following upper bound 

inf lim sup — log P„ ,(L n &B(u,5)) < — P(u). 



x,y)fj,(dx)fj,(dy). 
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Turning to the lower bound for (P n o L" 1 ), we show that for any \i G Ati(S) 

inf liminf — \ogP n (L n G B(u, 5)) > — k 2 [ [ log\x' y — y' y \n(dx)fj,(dy) + K [ logwdu. (14) 
<5>0 n->oo n 2 2 J J J 

Claim 3.6. We may assume, without loss of generality, that 

(i) n has no atoms 

(ii) S := supp(/i) is a compact subset of £ such that S D (J\f(w) U {0}) = 0. 
Proof. With the notations of the proof of the upper bound, we know that 

^n 2 J J log |x 7 — y' y \fi(dx)n(dy) + k J log w(x)fi(dx) = — J Fdfx® 2 

with F bounded from below. If fi has an atom, then J Fd/u,® 2 is infinite, and there is nothing 
to prove. Set 



4:= [-M]nsn |J }x--,x + -[\ 



and /i k := -^ii\ Ak . Then 

f Fdfi® 2 = fim / Fdn\ 



k— >oo 

Hence it suffices to prove ffl4"|) for fi k in the place of /i. Consequently, we may assume the 
support of \i to be contained in a finite union of compact intervals not meeting J\f(w) U {0}. 
This implies (ii). □ 

For j = 1, . . . ,p(n) let £ 3 - = £ ? - n ' ) be the P ^^~ J ' quantile of /i. Let 7 G N be as in above. 
Write £ 7 = (Q n) , . . . , ^7). Set £ p („)+i := inf S and £ := £1 + 1- Then, by Claim ESI 

-OO < Cp(n)+1 < £p(n) < • • • < £0 < OO. 

For 5 > 0, t G M p(n) write 

7T„(t) := {i = 1, . . . ,p{n) : tj > 0}, u n (t) := {1, . . . ,p(n)} \ 7r„(i), 

.(n) I 

i+il 



!„(*) :={< = !,..., p(n): |e?° - < *}, 



• J] n) (i, 5) := " h + 6] n E, j = 1, . . . ,p( 

7 W ff r\ ._ f fe, *j + 5] for j G 7r n (t), 
Jj l ' h 1 fo-<M;] for !/„(*), 

• Ut,5) ■.= Ylf:I^ n \t,s), Ut,sy : = Yi P M^ n \t,5y, 

Here, for M C E, we write M 7 := {m 7 : m G M}. 

Fix 5 > 0, and write G := 25). It follows from 0, Lemma 3.3] that for n large enough 
one has 

I„(£ (n) , 6) C {a; G £ p(n) : L n (x) G C7}. (15) 



n 
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Set y^ n) := (pf® := mi{w n (x) : x G [£,■ - 5, ^ + 5] U fj-i]}, j = 1, . . . ,p(n), and write 
ip n for the step function which equals tp^ on £j] and is zero elsewhere. 



For n as in (fT5l) we have 



= ~^{n) I Yl\ x i- X i\ fi X\w n (x] h ) n \xi\-^rm p{n) (dx). 

Observing that I n (£,5) 7 D J„(£V 7 ) and 

inf jn^O^T : x G I„(e,5) 7 | = mf jn™n(^) n : x G UMj > U/if^T 
we obtain that 

Pn(L n G G) > -L. JJ ((^(»)) n (| 6 | + 5)1-7) / JJ | X . _ Xj f mp(n) (^). 

Now, we can bound this last integral from below by 

J ( [0,5"*] « 7 ) x [- 5^ ,0] #"« (f 7 ) ) nRP(") •+ 



'([0,5^]#^(« 7 )x[-<5^,0]#"™(f 7 ))nIRP(™)'- ,,, , 
p(n)-l 



i<j— 1 i=l 
« p(n)-l 

x / i f |xj - 2 mp(n)(dx). (16) 

^([o,5^]#-«(e 7 )x[-<57,o]#^(« 7 ))nRp(»).+ ~ = "7 

By the change of variables u p c n ) = x p ( n )i u i~i — x % ~ x i-\ {i — p( n )i ■ ■ ■ , 2) one can bound 
ffTB"]) from below by 

P(») g ^ 2 ^p(n)-l ^ J7 ^ 5±3(p(n)-l)+l 

JJu? m p(n) (rfM) 



pRJ i=2 



So far, it has been shown that 
Pn{L n G G) 



p(n)-l 

if 



> n n 

i<j— 1 i=l 

1 // xn n / 9 \ PW" 1 / A7 \ ^ 2 (p( n )~ 1 ) + 1 



X 
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So we obtain 



1 lo g P n (L n eG) 



n 2 



>- (^)^EW-^l (17) 



+ (^^eW-ci , 18) 

+ ij>gp<"> (19) 

i 

+ ^E( 1 -^ 1 °g(^i+ 5 ) ( 20 ) 

i 

+ f-p(«)log7 + (p(n) - l)log-^-^J (21) 

+ ^( ^ + 2)( f dos(^) - lo S (p(n))). (22) 



Now, for (ITTj) and (TT8]) observe that 



log(?/ 7 - x 1 )fj,(dx)fj,(dy) 



x<y 



= E 

Kj 



log(?/ 7 - x 1 )fx(dx)fi(dy) 



(x,v)e[€ j+1 ,tj]x[Zi+i,ei] 



2 ~~~ J (x,y)e[£i+i&] x2 

£ ^ E - <w + ^ E Mff - ao. 



log|y 7 - x 7 \ji(dx)fx(dy) 
1 

2p(n) 



(12"U|) . (12"T|) . fT2"2l are easily seen to converge to zero. 

As to ( |T9l) . observe that ^ X]j=i logV 9 ^ = / log ^>„ef/i. Denote by Z a Lipschitz constant 
of log w on iS. For 77 > write 1% := max{| logw n (:r) — logw n (y)| : \x — y\ < r/} and define I 71 
analogously By (al), then, Q -> Z*. Note that Z" < Z77. Write M(n, 5) := \J je i n (S) [fj+i»fj n) ] 
and C := max{logw(x) : x G 5} — min{logw(x) : x G S}. 

Let e > 0. Since for all n and j > 1 one has fi ^j n ^]J = l/M 72 )) an d since 5 is fixed, 

one has fi(M(n, 5) c ) < e for large n. Now let n be large enough such that one also has 
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^25 _ p| < £ an( j || log^ _ log it; 1 1 < e. Then 



]2(S 



log1p n - log W\dfi < 

< 



j6/„(5) 



logV'n - logw n |d/i + e 

I logV'n - \ogw n \dfi + (C + 2e + l)e 



< p(n) — - 2/5 + (C + 2e + 3)e. 



So we have 



This implies ( !T4l) . 



lim sup 



71 .7=1 ^ 



0(5). 



Summing up, for any ji G .Mi(E) we have obtained 



inf lim sup — log P n (L n e B(fi,5)) < -H(fi) 



<5>0 



and 



inf lim inf — log P n {L ri e B(fx,5)) > -H(fj). 

8>0 n->oo n 

Using the exponential tightness of (P n o -L" 1 )^, we can apply [El Thm. 4.1.11] to obtain an 
LDP for (P n o L~ l ) n with rate H and speed n 2 . Setting A = G = A4i(S) in the lower and 
upper bound, we obtain 



lim — log Z n — - inf 

n-s-oo n A /LteMi(E) 

By Lemma [3.51 (ii), the right-hand side is < +oo. Now, 



1 1 

— logQ„(L n G A) = — (logP n (L n e A) - logZ r , 



for any Borel set A in A^i(S). Hence Theorem 13.11 is proven. 

4. Application to Gaussian Hamiltonian ensembles 



4.1. Gaussian Hamiltonian ensembles. Let 



j(n) c C d(n)xd(n) be the gpace of gQod 

Hamiltonians of symmetry class C. We wish to define Ad^- resp. Ad^-invariant probability 
measures on M.^ with the additional property that all matrix entries should be Gaussian 
and as many entries as possible should be independent. Plainly, what we are looking for is 
an analog for class C of the GOE, which was reviewed in the introduction and corresponds 
to C = AI. 

Let us look at an example: For C = CI, represent an element of as 

A 1 B -A 
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where A and B are real symmetric n x n matrices. The upper diagonal entries of A and B 
can be chosen as independent Gaussians. Then the distribution of X has a Lebesgue density 
of the form 

const x JJexp ( "^M II ex P ( ) ■ 

Now, a suitable choice of <7 a; ^, (Jb,k,i guarantees the invariance of the distribution of X. In 
fact, X being symmetric, 



Tr(X 2 ) = Tr(XX') = 2 Tr(AA' + BB') = 2 a% + foi 



2 



i<j i i<j 



hence 



8(T 2 v / ^ 2a 2 ^ 2(2<r 2 ) ^ 2a 2 ^ 2(2a 2 ) 

i<j i i<j i 

This means that in order to obtain an invariant distribution, we choose a 2 > and set 
a l,k,i = a b,k,i = (T 2 (k < I), o a ^ k = a b>k ,k = 2a 2 . We write GE^cr 2 ) for the probability 
distribution on obtained in this way. 

Note that the eigenvalues of X come in pairs ±A. So, if Ai, . . . , A n are the positive eigenvalues 
of X (under GE^((T 2 ), X has n distinct positive eigenvalues a.s.), then we have 

n 



8a 2 v ' 8a 2 ^ 1 Aa 2 
For comparison with the other symmetry classes, we write ipci = 8 and ipci = 4. 

(n) 

For a general symmetry class C, X e , we may proceed along the same lines. If the 
nonzero real and imaginary parts of the strictly upper triangular entries of (skew) symmetric 
or hermitian blocks have variance a 2 , then the nonzero real and imaginary parts of the 
diagonal entries of these blocks have variance 2<r 2 . This procedure determines an integer ip c 
such that the Lebesgue density of X has the form 

const, x exp f 3— Tr(X 2 ) 

V ^ccr 2 

If Ai,...,A p ( n ) are the positive eigenvalues without multiplicity of X (or all eigenvalues 
without multiplicity for C — A, AI, All), then there is an integer tp c defined by 

p(n) 

2 



y?co" 2 ip c o- 2 

Proceeding analogously to the example above and keeping in mind that X has 2-dimensional 
eigenspaces for C = All, CII and Dili, one obtains the following data: 
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c 


A 


AI 


All 


AIII 


B/D 


BDI 


Dili 


C 


CI 


CII 




4 


4 


8 


4 


4 


4 


8 


8 


8 


8 


ipc 


4 


4 


4 


2 


2 


2 


2 


4 


4 


2 



In the following subsections we apply our main theorem to the ensembles GE^ 



n ■ 



4.2. LDP for Wigner-Dyson ensembles A, AI, All. From Corollary 12.31 we see that 

one can subsume the joint eigenvalue density induced by GE^ fir)' C — A, AI, All, under 

the general form of (jHJ) by choosing p(n) = n, hence k = 1, 7 = 1, (3 = 1, 2, 4 according to 
C = AI, A, All, and 

1 



w n (xj) = exp 



'4a 2 Xj I ' 



(23) 



independent of n. Then, under GE^ f^J; {L n )n satisfies an LDP with good rate function 




log|x — y\ 1 jj(dx) fx(dy) + 



4a 2 



x 2 fx(dx) — const. 



(24) 



To apply the theory of logarithmic potentials with external fields, as presented in [24J, we 
exploit the fact that the unique minimizer /x* of Iq is also the unique minimizer of 

1 



Jc(jj) 




\og\x — y\ n(dx)n(dy) + 2 



Aa 2 (3 



x fj,(dx) 



One reads off from p. 284 (cf. p. 26 for notations) of [24J that fi* has Lebesgue density 



x- 



the density of Wigner's semicircle distribution of radius 2^Jo~ 2 f3. 

4.3. LDP for chiral ensembles BDI, AIII, CII. We will freely use the notation intro- 
duced above and in Section [21 We have ipc — 2 and p(n) = s(n) A t{n). For simplicity, we 
will assume that s(n) < t(n) for all iieN, hence a(n) = (3(t(n) — s(n)) + f3 — 1. Then we 

can subsume the joint eigenvalue density induced by GE^ f^r) under the general form of 

(12]) by setting 7 = 2 and 



w n (x 



/3(t(n)-s(n)) + /3-l _ x 2 

x « e 2^ 



hence 



w(x) = x^- 2 ^e-^. 



Then, by Theorem 13.11 the rate function of the LDP for Q n o L n is 




log 



x 2 — y 2 



■n(dx)n(dy) — k I log ( x 



,/3(1-2k) 



e 



fi(dx) — c, (25) 



for G A^i([0, oo[). For r > 0,a; > write T r (x) := x r . Now, /i* is the unique minimizer of 
Ic if, and only if, T 2 (/z*) is the unique minimizer of 
.2 







K 




log 



f - v\ 



Tv(dx)v(dy) — k I log (a/' 2 'e 2^ 



i/(cix) — c, 
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hence of 

J c {y):= J J log | 1 v(dx)v(dy) - 2 J log (x^e~^^) v(dx). (26) 

We can use the following facts from logarithmic potential theory with Laguerre weights: 
Lemma 4.1. For s > 0, A > the integral 

I.,xiy) ■= J J \o gl -^—p{dx)v{d y )-2 J log(x s e- Xx )u(dx), v G A4i([0, oo[) 
has a unique minimizer v* A with Lebesgue density 

l[a,b](x) ^y/(x-a)(b-x), (27) 

where 

a = a SjA = -(s + l- y/2s + l), b = 6 S>A = ~(s + 1 + y/2T+l). (28) 
A A 

Consequently, Ti(z/* A ) has Lebesgue density 

h^ x ) — V(^-a)(6-x2). (29) 
Proof. [21], IV (1.31), IV (5.18) □ 
We have 5=^ — 1 and A = 2g .lg K . Note that 2k < 1 by our assumptions. This yields 

a = 2a 2 p Q - v/ K (l - , b = 2a 2 (3 Q + - k)) , (30) 

and from fl29|) we conclude that the minimizer fi* of Ic has Lebesgue density 

i^^) ^L^ V{x2 ~ a){b ~ x2) (31) 

with a, 6 as in (1301). 



Remark 4.2. T2(/x*) is related to, but does not coincide with, the Marcenko-Pastur dis- 
tribution, which arises as limit of the empirical eigenvalue measure of Wishart (or more 
generally, sample covariance) matrices. Details on this relationship can be found in the com- 
panion paper [2TJ. The LDP for Wishart matrices, first proven by Hiai and Petz in [Tj5], is 
a consequence of Theorem 13.11 above. 

4.4. LDP for BdG ensembles B, D, Dili, C, CI. Although there are five BdG ensem- 
bles, there are only four natural large n limits for these ensembles. This is because B and 
D are iso(n) for odd resp. even n, and, as n tends to infinity, these series should be thought 
of as interlaced. It is consistent to do so, because we will see that the parameter a, where 
B and D differ (see table in Corollary 12.31) . does not affect the rate function. By the same 
token, it is in fact consistent to consider one rather than two large n limits for Dili, because 
the even and odd cases of Dili only differ w.r.t. a. 

We subsume the joint eigenvalue density induced by GE^ 1 -' (^j, C = B, D, Dili, C, CI 
under the general form (j3J) by setting 7 = 2, (a, 0) according to the table after Corollary 
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12.31 p(n) = [f J fc> r C = Dili and p(ri) = n otherwise, hence k — | or k — 1, respectively. 
This yields 

w n (x) = x°/ ra exp( — p— - x 2 ^ 

hence 



1 



wCx) = exp( x ) 

with ^b/d = 2 and = V'ci = V'diii = 4. 

Then, under GE^ (^~)' (-^n)n satisfies an LDP with good rate function 



IM = t^ 2 J J log |a: 2 — j/ 2 1 1 fi(dx) fi(dy) + — - 2 ^ x 2 /j(dx) — const. 



To determine the unique minimizer /i* of Iq, one can proceed as in Subsection 14.31 and apply 
Lemma I4TT1 with s = and A = ^ ^2^ K - By /x* is a quarter circle distribution, given by 
the Lebesgue density 

2 
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